Genomic Disorders in CKD across the Lifespan.

SIGNIFICANCE STATEMENT: Pathogenic structural genetic variants, also known as genomic disorders, have been associated with pediatric CKD. This study extends those results across the lifespan, with genomic disorders enriched in both pediatric and adult patients compared with controls. In the Chronic Renal Insufficiency Cohort study, genomic disorders were also associated with lower serum Mg, lower educational performance, and a higher risk of death. A phenome-wide association study confirmed the link between kidney disease and genomic disorders in an unbiased way. Systematic detection of genomic disorders can provide a molecular diagnosis and refine prediction of risk and prognosis. BACKGROUND: Genomic disorders (GDs) are associated with many comorbid outcomes, including CKD. Identification of GDs has diagnostic utility. METHODS: We examined the prevalence of GDs among participants in the Chronic Kidney Disease in Children (CKiD) cohort II ( n =248), Chronic Renal Insufficiency Cohort (CRIC) study ( n =3375), Columbia University CKD Biobank (CU-CKD; n =1986), and the Family Investigation of Nephropathy and Diabetes (FIND; n =1318) compared with 30,746 controls. We also performed a phenome-wide association analysis (PheWAS) of GDs in the electronic MEdical Records and GEnomics (eMERGE; n =11,146) cohort. RESULTS: We found nine out of 248 (3.6%) CKiD II participants carried a GD, replicating prior findings in pediatric CKD. We also identified GDs in 72 out of 6679 (1.1%) adult patients with CKD in the CRIC, CU-CKD, and FIND cohorts, compared with 199 out of 30,746 (0.65%) GDs in controls (OR, 1.7; 95% CI, 1.3 to 2.2). Among adults with CKD, we found recurrent GDs at the 1q21.1, 16p11.2, 17q12, and 22q11.2 loci. The 17q12 GD (diagnostic of renal cyst and diabetes syndrome) was most frequent, present in 1:252 patients with CKD and diabetes. In the PheWAS, dialysis and neuropsychiatric phenotypes were the top associations with GDs. In CRIC participants, GDs were associated with lower serum magnesium, lower educational achievement, and higher mortality risk. CONCLUSION: Undiagnosed GDs are detected both in children and adults with CKD. Identification of GDs in these patients can enable a precise genetic diagnosis, inform prognosis, and help stratify risk in clinical studies. GDs could also provide a molecular explanation for nephropathy and comorbidities, such as poorer neurocognition for a subset of patients.

Genetic regulation of serum IgA levels and susceptibility to common immune, infectious, kidney, and cardio-metabolic traits.

Immunoglobulin A (IgA) mediates mucosal responses to food antigens and the intestinal microbiome and is involved in susceptibility to mucosal pathogens, celiac disease, inflammatory bowel disease, and IgA nephropathy. We performed a genome-wide association study of serum IgA levels in 41,263 individuals of diverse ancestries and identified 20 genome-wide significant loci, including 9 known and 11 novel loci. Co-localization analyses with expression QTLs prioritized candidate genes for 14 of 20 significant loci. Most loci encoded genes that produced immune defects and IgA abnormalities when genetically manipulated in mice. We also observed positive genetic correlations of serum IgA levels with IgA nephropathy, type 2 diabetes, and body mass index, and negative correlations with celiac disease, inflammatory bowel disease, and several infections. Mendelian randomization supported elevated serum IgA as a causal factor in IgA nephropathy. African ancestry was consistently associated with higher serum IgA levels and greater frequency of IgA-increasing alleles compared to other ancestries. Our findings provide novel insights into the genetic regulation of IgA levels and its potential role in human disease.

Copy Number Variant Analysis and Genome-wide Association Study Identify Loci with Large Effect for Vesicoureteral Reflux.

BACKGROUND: Vesicoureteral reflux (VUR) is a common, familial genitourinary disorder, and a major cause of pediatric urinary tract infection (UTI) and kidney failure. The genetic basis of VUR is not well understood. METHODS: A diagnostic analysis sought rare, pathogenic copy number variant (CNV) disorders among 1737 patients with VUR. A GWAS was performed in 1395 patients and 5366 controls, of European ancestry. RESULTS: Altogether, 3% of VUR patients harbored an undiagnosed rare CNV disorder, such as the 1q21.1, 16p11.2, 22q11.21, and triple X syndromes ((OR, 3.12; 95% CI, 2.10 to 4.54; P=6.35×10-8) The GWAS identified three study-wide significant and five suggestive loci with large effects (ORs, 1.41-6.9), containing canonical developmental genes expressed in the developing urinary tract (WDPCP, OTX1, BMP5, VANGL1, and WNT5A). In particular, 3.3% of VUR patients were homozygous for an intronic variant in WDPCP (rs13013890; OR, 3.65; 95% CI, 2.39 to 5.56; P=1.86×10-9). This locus was associated with multiple genitourinary phenotypes in the UK Biobank and eMERGE studies. Analysis of Wnt5a mutant mice confirmed the role of Wnt5a signaling in bladder and ureteric morphogenesis. CONCLUSIONS: These data demonstrate the genetic heterogeneity of VUR. Altogether, 6% of patients with VUR harbored a rare CNV or a common variant genotype conferring an OR >3. Identification of these genetic risk factors has multiple implications for clinical care and for analysis of outcomes in VUR.

The genetic architecture of membranous nephropathy and its potential to improve non-invasive diagnosis.

Membranous Nephropathy (MN) is a rare autoimmune cause of kidney failure. Here we report a genome-wide association study (GWAS) for primary MN in 3,782 cases and 9,038 controls of East Asian and European ancestries. We discover two previously unreported loci, NFKB1 (rs230540, OR = 1.25, P = 3.4 × 10-12) and IRF4 (rs9405192, OR = 1.29, P = 1.4 × 10-14), fine-map the PLA2R1 locus (rs17831251, OR = 2.25, P = 4.7 × 10-103) and report ancestry-specific effects of three classical HLA alleles: DRB1*1501 in East Asians (OR = 3.81, P = 2.0 × 10-49), DQA1*0501 in Europeans (OR = 2.88, P = 5.7 × 10-93), and DRB1*0301 in both ethnicities (OR = 3.50, P = 9.2 × 10-23 and OR = 3.39, P = 5.2 × 10-82, respectively). GWAS loci explain 32% of disease risk in East Asians and 25% in Europeans, and correctly re-classify 20-37% of the cases in validation cohorts that are antibody-negative by the serum anti-PLA2R ELISA diagnostic test. Our findings highlight an unusual genetic architecture of MN, with four loci and their interactions accounting for nearly one-third of the disease risk.

Author Correction: The copy number variation landscape of congenital anomalies of the kidney and urinary tract.

In the version of this article initially published, affiliation 38 incorrectly read "ICNU-Nephrology and Urology Department, Barcelona, Spain"; "Renal Division, Hospital Clinic, IDIBAPS, University of Barcelona, Barcelona, Spain" is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

The Burden of Candidate Pathogenic Variants for Kidney and Genitourinary Disorders Emerging From Exome Sequencing.

Background: Exome sequencing is increasingly being used for clinical diagnostics, with an impetus to expand reporting of incidental findings across a wide range of disorders. Analysis of population cohorts can help reduce risk for genetic variant misclassification and resultant unnecessary referrals to subspecialists. Objective: To examine the burden of candidate pathogenic variants for kidney and genitourinary disorders emerging from exome sequencing. Design: Secondary analysis of genetic data. Setting: A tertiary care academic medical center. Patients: A convenience sample of exome sequence data from 7974 self-declared healthy adults. Measurements: Assessment of the prevalence of candidate pathogenic variants in 625 genes associated with Mendelian kidney and genitourinary disorders. Results: Of all participants, 23.3% carried a candidate pathogenic variant, most of which were attributable to previously reported variants that had implausibly high allele frequencies. In particular, 25 genes (discovered before the creation of the Exome Aggregation Consortium, a genetic database comprising data from a large control population) accounted for 67.7% of persons with candidate pathogenic variants. After stringent filtering based on allele frequency, 1.4% of persons still had a candidate pathogenic variant, an excessive rate given the prevalence of monogenic kidney and genitourinary disorders. Manual annotation of a subset of variants showed that the majority would be classified as nonbenign under current guidelines for clinical sequence interpretation and could prompt subspecialty referrals if returned. Limitation: Limited access to health record data prevented comprehensive assessment of the phenotypic concordance with genetic diagnoses. Conclusion: Widespread reporting of incidental genetic findings related to kidney and genitourinary disorders will require stringent curation of clinical variant databases and detailed case-level review to avoid genetic misdiagnosis and unnecessary referrals. These findings motivate similar analyses for genes relevant to other medical subspecialties. Primary Funding Source: National Institute of Diabetes and Digestive and Kidney Diseases and National Human Genome Research Institute.

The copy number variation landscape of congenital anomalies of the kidney and urinary tract.

Congenital anomalies of the kidney and urinary tract (CAKUT) are a major cause of pediatric kidney failure. We performed a genome-wide analysis of copy number variants (CNVs) in 2,824 cases and 21,498 controls. Affected individuals carried a significant burden of rare exonic (that is, affecting coding regions) CNVs and were enriched for known genomic disorders (GD). Kidney anomaly (KA) cases were most enriched for exonic CNVs, encompassing GD-CNVs and novel deletions; obstructive uropathy (OU) had a lower CNV burden and an intermediate prevalence of GD-CNVs; and vesicoureteral reflux (VUR) had the fewest GD-CNVs but was enriched for novel exonic CNVs, particularly duplications. Six loci (1q21, 4p16.1-p16.3, 16p11.2, 16p13.11, 17q12 and 22q11.2) accounted for 65% of patients with GD-CNVs. Deletions at 17q12, 4p16.1-p16.3 and 22q11.2 were specific for KA; the 16p11.2 locus showed extensive pleiotropy. Using a multidisciplinary approach, we identified TBX6 as a driver for the CAKUT subphenotypes in the 16p11.2 microdeletion syndrome.

Diagnostic Utility of Exome Sequencing for Kidney Disease.

BACKGROUND: Exome sequencing is emerging as a first-line diagnostic method in some clinical disciplines, but its usefulness has yet to be examined for most constitutional disorders in adults, including chronic kidney disease, which affects more than 1 in 10 persons globally. METHODS: We conducted exome sequencing and diagnostic analysis in two cohorts totaling 3315 patients with chronic kidney disease. We assessed the diagnostic yield and, among the patients for whom detailed clinical data were available, the clinical implications of diagnostic and other medically relevant findings. RESULTS: In all, 3037 patients (91.6%) were over 21 years of age, and 1179 (35.6%) were of self-identified non-European ancestry. We detected diagnostic variants in 307 of the 3315 patients (9.3%), encompassing 66 different monogenic disorders. Of the disorders detected, 39 (59%) were found in only a single patient. Diagnostic variants were detected across all clinically defined categories, including congenital or cystic renal disease (127 of 531 patients [23.9%]) and nephropathy of unknown origin (48 of 281 patients [17.1%]). Of the 2187 patients assessed, 34 (1.6%) had genetic findings for medically actionable disorders that, although unrelated to their nephropathy, would also lead to subspecialty referral and inform renal management. CONCLUSIONS: Exome sequencing in a combined cohort of more than 3000 patients with chronic kidney disease yielded a genetic diagnosis in just under 10% of cases. (Funded by the National Institutes of Health and others.).

Exome-wide Association Study Identifies GREB1L Mutations in Congenital Kidney Malformations.

Renal agenesis and hypodysplasia (RHD) are major causes of pediatric chronic kidney disease and are highly genetically heterogeneous. We conducted whole-exome sequencing in 202 case subjects with RHD and identified diagnostic mutations in genes known to be associated with RHD in 7/202 case subjects. In an additional affected individual with RHD and a congenital heart defect, we found a homozygous loss-of-function (LOF) variant in SLIT3, recapitulating phenotypes reported with Slit3 inactivation in the mouse. To identify genes associated with RHD, we performed an exome-wide association study with 195 unresolved case subjects and 6,905 control subjects. The top signal resided in GREB1L, a gene implicated previously in Hoxb1 and Shha signaling in zebrafish. The significance of the association, which was p = 2.0 × 10-5 for novel LOF, increased to p = 4.1 × 10-6 for LOF and deleterious missense variants combined, and augmented further after accounting for segregation and de novo inheritance of rare variants (joint p = 2.3 × 10-7). Finally, CRISPR/Cas9 disruption or knockdown of greb1l in zebrafish caused specific pronephric defects, which were rescued by wild-type human GREB1L mRNA, but not mRNA containing alleles identified in case subjects. Together, our study provides insight into the genetic landscape of kidney malformations in humans, presents multiple candidates, and identifies SLIT3 and GREB1L as genes implicated in the pathogenesis of RHD.

Sample preparation method considerations for integrated transcriptomic and proteomic analysis of tumors.

Sample processing protocols that enable compatible recovery of differentially expressed transcripts and proteins are necessary for integration of the multiomics data applied in the analysis of tumors. In this pilot study, we compared two different isolation methods for extracting RNA and protein from laryngopharyngeal tumor tissues and the corresponding adjacent normal sections. In Method 1, RNA and protein were isolated from a single tissue section sequentially and in Method 2, the extraction was carried out using two different sections and two independent and parallel protocols for RNA and protein. RNA and protein from both methods were subjected to RNA-seq and iTRAQ-based LC-MS/MS analysis, respectively. Analysis of data revealed that a higher number of differentially expressed transcripts and proteins were concordant in their regulation trends in Method 1 as compared to Method 2. Cross-method comparison of concordant entities revealed that RNA and protein extraction from the same tissue section (Method 1) recovered more concordant entities that are missed in the other extraction method (Method 2) indicating heterogeneity in distribution of these entities in different tissue sections. Method 1 could thus be the method of choice for integrated analysis of transcriptome and proteome data.